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Abstract. Feedback can be utilized to convert information into useful work, making 
it an effective tool for increasing the performance of thermodynamic engines. Using 
feedback reversibility as a guiding principle, we devise a method for designing optimal 
feedback protocols for thermodynamic engines that extract all the information gained 
during feedback as work. Our method is based on the observation that in a feedback- 
reversible process the measurement and the time-reversal of the ensuing protocol both 
prepare the system in the same probabilistic state. We illustrate the utility of our 
method with two examples of the multi-particle Szilard engine. 
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1. Introduction 

An important application of feedback is to increase the performance of thermodynamic 
engines by converting the information gathered during feedback into mechanical 
work [H El m HI [5l ini [7]. However, for feedback implemented discretely - through 
a series of feedback loops initiated at predetermined times - the second law of 
thermodynamics for discrete feedback limits the maximum amount of work that can 
be extracted [HI El HDl [HI [121 [iSl [E]. Namely, the average work extracted (W) during 
a thermodynamic process with discrete feedback in which a system is driven from one 
equilibrium state at temperature T to another equilibrium state at the same temperature 
is bounded by the difference between the information gained during feedback (/) and 
the average free energy difference (AF): 

{W)<kT{I)-{AF), (1) 

where k is Boltzmann's constant. Here, (/) is the mutual information between the 
microscopic state of the system and the measurement outcomes, and (AF) is the average 
free energy difference between the initial equilibrium state and the final equilibrium 
state, which may differ for each measurement outcome. Notice is expressed in terms 
of the extracted work, since we have in mind applications to thermodynamic engines. 
This differs from the more common convention of using the work done on the system, 
which is minus the work extracted [8l [9l [Tnl [TTl [T2l [T3l [Tl] . 

Optimal thermodynamic engines extract the maximum amount of work, saturating 
the bound in ([1]) [{W) = kT{I) — (AF)]. Their design often proceeds in two steps. 
One first selects a physical observable M to be measured. Then, associated to each 
measurement outcome m, one chooses a unique protocol for varying a set of external 
parameters A during a time interval from t = to r. A'" = {AJ"}[^o- For the process to 
be optimal the collection of protocols {A"*} must be designed to extract as work all the 
information gained from the measurement. 

While at first it may not be obvious how to design a collection of optimal 
protocols O [6], there is a generic procedure for constructing such a collection given 
a physical observable M [HI [IHl [161 [13 [13 [IH] ; specifically, the optimal protocol is to 
instantaneously switch the Hamiltonian immediately after the measurement - through 
an instantaneous change of the external parameters - so that the probabilistic state 
of the system conditioned on the measurement outcome is an equilibrium Boltzmann 
distribution with respect to the new Hamiltonian. The external parameters are then 
reversibly adjusted to their final value, completing the protocol. While such a protocol 
can always be constructed theoretically, it may be difficult to realize experimentally: 
one may need access to an infinite number of external parameters in order to affect 
the instantaneous switching of the Hamiltonian [19]. Furthermore, there are optimal 
protocols that cannot be constructed by implementing this generic procedure. Hence, 
it is worthwhile to develop alternative procedures for engineering collections of optimal 
protocols. 
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In a recent article, we characterized optimal feedback processes, demonstrating that 
they are feedback reversible - indistinguishable from their time- reversals [20]. There we 
pointed to the possibility of exploiting feedback reversibility in the design of optimal 
thermodynamic engines. In this article, we take the next step by explicitly formulating a 
recipe for engineering a collection of optimal feedback protocols for a given observable M 
using feedback reversibility as a guiding principle. We present our method in section |2l 
generalizing the generic procedure outlined in the previous paragraph. We then illustrate 
our method in section |3] with two pedagogical models inspired by the multi-particle 
Szilard engine recently introduced in [5], and subsequently analyzed in [21]: a classical 
two-particle Szilard engine with hard-core interactions, and a classical iV-particle Szilard 
engine with short-ranged, repulsive interactions. In each model, we design a different 
collection of feedback protocols, demonstrating the utility and versatility of our method. 
Concluding remarks are offered in section H] with a view towards potential applications 
of our method to quantum feedback. 

2. Measurement and preparation 

In this section, we describe a general method for designing optimal feedback protocols. 
Our analysis is based on a theoretical framework characterizing the thermodynamics of 
feedback formulated in[9l[T0l[IIl[T2l[T3l[Tll[20]. 

Consider a classical system whose position in phase space at time t is Zt- The 
system, initially in equilibrium at temperature T, is driven by varying a set of external 
control parameters A initially at Aq from time t = to r using feedback. At time 
t = tm, a.n observable M is measured whose outcomes m occur randomly with probability 
P{m\zt^) depending only on the state of the system at the time of measurement zt^. The 
protocol, denoted as A™ = {A™}[^q, depends on the measurement outcome after time 
tm- Thermal fluctuations cause the system to trace out a random trajectory through 
phase space 7 = {zt}l=Q. The work extracted along this trajectory is W['j; A'"], and the 
reduction in our uncertainty due to the measurement is [H [101 120] 



where P{m) is the probability of obtaining measurement outcome m. For error-free 
measurements, which we consider in our illustrative examples below, the measurement 
outcome is uniquely determined by the state of the system at the time of measurement. 
Consequently, P{m\zt^) is always either zero or one. When P{m\zt^) = 1, ([2]) reduces 
to 



When P{m\zt^) = 0, (|2]) is divergent; however, this divergence occurs with zero 
probability, and therefore does not contribute to the average in ([T]). Finally, the change 
in free energy from the initial equilibrium state, -F(Ao), to the final equilibrium state, 




(2) 




(3) 
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F(A™), denoted as AF[A™] = F(A^) — -F(Ao), is realization dependent, since the final 
external parameter value at time r depends on the measurement outcome m. 

Associated to the feedback process is a distinct thermodynamic process called the 
reverse process [TOl [131 ISO]- The reverse process begins by first randomly selecting a 
protocol A"^ according to P{m). The system is then prepared in an equilibrium state 
at temperature T with external parameters set to A™. From time t = to r, the 
system is driven by varying the external parameters according to the time-reversed 
conjugate protocol A"* = {\t}l=Q, where AJ" = A^^. For every trajectory 7 = {zt}l^Q 
of the forward process there is a time-reversed conjugate trajectory 7 = {5t}[^Q, where 
Zt = z*_f- and * denotes momentum reversal. 

A feedback process that is indistinguishable from its reverse process is called 
feedback reversible [20]. A useful microscopic expression for the present considerations is 
in terms of the phase space densities along the feedback process and the corresponding 
reverse process. Namely, the phase space density of the feedback process at time t 
conditioned on executing protocol A™, pi^ztlA"^), is identical to the phase space density 
in the reverse process at time t — t conditioned on executing protocol A*", p{zr-t\-A"^)'- 

p{zt\An = pizr-^ti^n- (4) 

Additionally, 

ly [7, A'"] = kTI[-f, A™] - AF[A"*] (5) 

for every realization [20]. For cyclic {AF = 0) feedback-reversible processes, such as our 
illustrative examples, ([5]) is simply iy[7. A™] = A;T/[7, A*"]. 

We now utilize (jl]) and ([5]) to develop a method for designing optimal feedback 
processes (or equivalently feedback-reversible processes). Our method is based on the 
observation that (jl]) has a noteworthy interpretation at the measurement time t = t^'- 

p{ztjAn = pr^r-tj~An- (6) 

Specifically, p{zt^\A"^) is the phase space density of the system at the time of the 
measurement conditioned on implementing protocol A™; it represents our knowledge 
about the microscopic state of the system immediately after the measurement. We 
therefore refer to it as the post-measurement state. The right hand side of ([6]), 
p(5t— |A™), is the phase space density at time t = t — tm produced by the reverse 
process when protocol A*" is executed; it is the probabilistic state of the system prepared 
(or produced) by using protocol A™ in the reverse process. Thus, we refer to p{zr-t^\A^) 
as the prepared state. With this terminology, ([6]) states that for a process to be 
feedback reversible the state prepared by the reverse process must be identical to the 
post-measurement state. This insight is our main tool for designing optimal feedback 
protocols. Instead of focusing on the feedback process, we search for a protocol that 
prepares the post-measurement state. We call this procedure preparation. Once we have 
chosen our protocols, we can verify their effectiveness by checking the equality in ([5]); 
the deviation from equality in ([5]) is a measure of the the reversibility of each of the 
protocols in {A"^}. 
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3. Applications to the multi-particle Szilard engine 

In this section, we apply the preparation method presented in section [2] to two classical 
extensions of the Szilard engine inspired by the quantum multi-particle Szilard engine 
considered by Kim et. al. in [5]. In section 13.11 we design a collection of optimal 
protocols for a classical Szilard engine composed of two square particles with hard- 
core interactions. An A-particle Szilard engine consisting of ideal point particles with 
short-ranged, repulsive interactions is analyzed in section 13.21 In both examples, we 
verify that our protocols are optimal through analytic calculations of the work and 
information. 

3.1. Two-particle Szilard engine 

To illustrate the utility of our method, we now analyze a two-particle Szilard engine. 
We have in mind two indistinguishable square hard-core particles with linear dimension 
d confined to a two-dimensional box of width and height Ly, pictured in figure [1] 
The particles have a hard-core interaction with the walls, entailing that the center of 




L 



Figure 1. Depiction of the two-partiele Szilard engine composed of two square hard- 
core particles of width d in a two-dimensional box of width and height Ly. Feedback 
protocols are initiated by infinitely slowly inserting a partition from the bottom edge 
of the box, as indicated by the vertical gray arrow, dividing the box into equal halves 
of width Lx/2. 

the particles must be at least a distance d/2 from the walls. The box is in weak thermal 
contact with a thermal reservoir at temperature kT = 1. 

Work is extracted using a cyclic, isothermal feedback protocol performed infinitely 
slowly, as illustrated in figure [2j Since the process is cyclic, (AF) = 0, and we only need 
to investigate the extracted work. In addition, since the process is infinitely slow and 
isothermal, the work can be expressed in terms of partition functions, as in [22]. There 
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are two configurational partition functions that will prove useful: the first, denoted 
Z2{x,y), is the partition function for the state when both particles are in the same 
box of width x and height y; the second, Z2{x,y), is the partition function for the 
state where the particles are in seperate boxes, each of width x and height y. The 
calculation of these partition functions is a straightforward though lengthy exercise in 



integral calculus, which we outline in Appendix A 



We initiate the feedback protocol with the engine in thermal equilibrium at 
temperature kT = 1. We then infinitely slowly insert a thin partition from below, 
dividing the box into two equal halves along the horizontal direction, as depicted in 
figure [H Because the particles are hard-bodied and of finite size, the insertion of 




Figure 2. Illustration of the three protocols executed in the two-particle Szilard engine 
associated to the three measurement outcomes A, B, and C. 

the partition extracts work. As we slowly insert the partition, the system remains 
in equilibrium and able to explore its entire phase space until the leading tip of the 
partition is one particle length d from the box's top wall. At which point, the particles 
are too large to pass between the left and right half of the box. At that moment, each 
particle becomes trapped in one half of the box; either they both become trapped in 
the same half of the box, or each is trapped in a separate half of the box. The partition 
function at that moment, being a sum over all distinct microscopic configurations, is 
then the sum of the partition function when they both become trapped in the left 
(or right) half, Z2{Lx/2, Ly), plus the partition function when they become trapped in 
separate halves, Z2(L^/2, Ly): 2Z2{Lx/2, Ly) + Z2(I/a./2, Ly). The work extracted up to 
that instant is determined from the ratio of the partition function at that moment to 
the initial partition function Z2{Lx,Ly) as 

Wps,rtiLa;,Ly) = In 



2Z2{Lj2,Ly) + Z2{Lj2,L, 

Z2{Lx, Ly) 



(7) 



Once the distance between the leading tip of the partition and the far wall of the box is 
less than d, neither particle is able to fit in the space between the tip and the wall. The 



Designing optimal feedback engines 



7 



partition's tip is no longer able to push on the particles, and as a result no additional 
work beyond that in ([7]) is extracted. 

Next, we measure in which half of the box the two particles are located. There are 
three outcomes, which we label A, B, and C, see figure [2l Outcomes A and C occur 
when both particles are found in the same half of the box, whereas outcome B occurs 
when each particle is found in a separate half of the box. Since the partition functions 
Z2 and Z2 count the number of distinct microscopic configurations, we can express the 
change in uncertainties associated to each outcome by inserting these partition functions 
into (Bl): 



Ia 



IB 



-In 



-In 



2Z2{Lj2,Ly) + Z2{Lj2,Ly)_ 
Z2{Lj2,Ly) 



(8) 

2Z2{Lj2,Ly) + Z2{Lj2,Ly)\ ' 

If both particles are found in the same half of the box (outcome A or C), the 
optimal protocol is to quasi-statically shift the partition to the opposite end of the box, 
as in the single-particle Szilard engine [23], extracting work 



shift 



In 



Z2{L. 



_^2(-^x/2, Ly) ' 

Summing ([7j) and (|TOl) . we find that the work extracted during the feedback protocol 
associated to measurement outcome A (or C) is 

Wa = W^,rt{Lx, Ly) + Wshift (11) 

^ 1^ \ 2Z2iLj2,Ly) + Z2{Lx/2,Ly) } ^ ^^2) 
Z2{Lx/2,Ly) J 

which equals in ([8]). Thus, according to ([5]) this protocol is optimal as expected, 
since this protocol when run in reverse clearly prepares the post-measurement state 
conditioned on A. 

When each particle is found in a separate half of the box (outcome B), the optimal 
protocol is less clear. The motion of the piston in either direction requires work rather 
than extracts it. Kim et. ai, for instance, opt to extract the partition without obtaining 
any useful work [5]: the information in the measurement is wasted. However, our 
discussion in section |2] suggests a way to design an optimal cyclic protocol: the protocol 
must drive the system from the state post measurement outcome B back to the initial 
state and when run in reverse must prepare the state associated to outcome B by 
segregating each particle into a different half of the box. When the particles do not 
interact, there is no obvious optimal protocol. However, in our model we can exploit the 
particle interactions. Specifically due to the hard-core interactions, there is a greater 
likelihood of trapping the particles in separate halves of the box upon inserting the 
partition when the box is smaller. This observation suggests the following protocol 
executed in response to measurement outcome B. 
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After the partition is inserted, we infinitely slowly compress the box until its width 
is Ix > 2d and its height is ly > d. The extracted work during compression is 

^2(^/2, ly) 



comp 



In 



Z2{Lj2,Ly) 

Next, the partition is removed infinitely slowly, extracting 



-W^part(/: 



Finally, the box is expanded back to its original size extracting 



exp 



In 



XI ly) 



(13) 
[see ([7])] work. 

(14) 



Combining the sum of (ED, (III), (111), and -IVpart(^x, ly), with (ED, 
algebraic manipulation, that the deviation from reversibility [cf. 

as 

, ^^2(/./2,/,) 



Wn-I 



B 



-In 



we find, after a simple 
can be expressed 

'Ulx/2,ly) ■ ^^^^ 

Note that Wb — Ib only depends on the size of the compressed box with dimensions 
Ix X ly. To investigate the reversibility of our protocol, we study the dependence of 
Wb — Ib on the compressed box size. To simplify our analysis, we only consider boxes 
such that Ix = 2ly. In figure |3l we plot Wb — /b as a function of the box size parameter 
^ = Ix/d = 2ly/d. The smaller ^ the smaller the box. Notice that Wb — Ib < 0- We also 




Figure 3. Plot of the deviation from reversibility Wb — Ib for the two- particle Szilard 
engine protocol implemented in response to measuring each particle in a separate half 
of the box (outcome B) as a function of the box size parameter ^ = Ix/d = 2ly/d. 

observe that the process becomes reversible {Wb — /b = 0) when ,^ < 4 (/^ < 4(i and 
ly < 2d)] the box is so small when < 4 that both particles cannot fit into the same 
half of the box. Consequently, when the partition is inserted during the reverse process 
each particle is confined to a separate half of the box, preparing the post-measurement 
state with probability one. 
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To confirm that our protocol can be optimal, we plot in figure S] the total average 
work extracted {W) = PaWa+PbWb+PcWc - where Pj is the probability to implement 
protocol j = A, B,C - as a. function of the box size parameter ^. Again, we see that 




3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 



Figure 4. Plot comparing the total average work extracted (W) (solid) in our two- 
particle Szilard engine protocol to the information (J) (long dashed) as a function 
of the box size parameter ^ = Ix/d = 2ly/d for a box initially of size ~ 20d by 
Ly = lOd. For reference, we have included the total average work that would have 
been extracted using the protocol introduced by Kim et al. in Ref. [5], (Wk) (dashed). 

when ^ < 4 our protocol becomes optimal: (W) = (/). For comparison, we have 
included in figure H] the work extracted when implementing the protocol proposed in 
Ref. [5], (Wk), where the partition is slowly removed in response to outcome B. 

Further insight can be gained by noting that the ratio Z2/ Z2 in f|T5|) . which controls 
the degree of reversibility, has a simple physical interpretation in terms of the change 
in free energy during an irreversible mixing of two indistinguishable particles, each in 
separate boxes of sizes 1^/2 x ly, into one box of the same size, 1^/2 x ly-. 

^ Z2{lj2,ly) 

_Z2{lx/2, ly) 

Thus, this protocol is reversible when there is an infinite free energy difference between 
the states in which both particles are in the same box and where each particle is 
in a separate box. For an ideal gas AF^i^ = In 2: two indistinguishable ideal gas 
particles confined to the same box have half as many distinct microscopic configurations 
than when they are in seperate boxes. For ideal gases our protocol is not optimal 
(AFmix 7^ 00 and Wb — 7^ 0), as it exploits particle interactions. Nevertheless, there 
may exist other protocols that are optimal for ideal gases. In particular, such a collection 
could be devised using the generic procedure outlined in the Introduction, where the 
Hamiltonain is instantaneously switched immediately after the measurement so that 
the post-measurement state is described by an equilibrium Boltzmann distribution with 



AFmix = - In 



(16) 
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respect to the new Hamiltonian [HI [151 [ISl [13 [ISl [IH] ; however, this new Hamiltonian 
would contain an interaction potential that forces the particles to segregate themselves 
into opposite halves of the box. 

3.2. N -particle Szilard engine 

As a final illustration, we present an optimal feedback protocol for a classical A^-particle 
Szilard engine. Consider indistinguishable, classical, point particles with short- 
ranged, repulsive interactions confined to a box of volume V in weak thermal contact 
with a thermal reservoir at temperature kT = 1. The protocol begins by quickly and 
isothermally inserting an infinitely thin partition into the box dividing it into two equal 
halves of volume V/2. Since this is performed rapidly and the particles are infinitely 
small, the particles never have an opportunity to interact with the partition implying 
that this insertion requires no work. We then measure the number of particles in the 
left half of the box. Based on the outcome, we implement a cyclic, isothermal feedback 
protocol. 

The change in uncertainty when n particles are found in the left half of the box 
(A^ — n particles in the right half) is, from (|3]), 

r 1 A^' 1 

4 = -In T^-TTV-r • (17) 

2^^ n!(A — nj! J 

This information can be extracted completely as work by implementing the following 
protocol. First, we slowly lower n [N — n) localized potential minima or trapping 
potentials to a depth E in the left (right) half of the box. The trapping potentials are 
assumed to be deep compared to the thermal energy {E ^ kT), but shallow compared to 
the interaction energy; so that only one particle is confined in each trapping potential, as 
depicted in figure [3 The partition is then quickly removed, and the trapping potentials 
are slowly turned off. 

Work is only extracted when the trapping potentials are turned on or off. Since 
these processes are very slow, the work extracted can be computed in terms of partition 
functions. Assuming that the volume V of the box is large compared with the interaction 
length, we can approximate the configurational partition function for the equilbrium 
state prior to inserting the partition as 

yN 

ZiV) = -. (18) 

After making the measurement and finding n particles in the left half of the box, the 
configurational partition function is 

^"'^)-S(]vb^(T)"' <^^' 

After lowering the trapping potentials to a depth E each particle is confined to a unique 
trapping potential of volume v. At which point the configurational partition function is 

Z^{v) = v^e-^^. (20) 
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Figure 5. Depiction of the 3- particle Szilard engine protocol in a box of volume V 
after having turned on the trapping potentials when 2 particles were found in the left 
half the box. Each of the three particles, pictured as black circles, is confined to a 
distinct trapping potential of volume w, pictured as a dashed circle. 



In terms of these partion functions, the work extracted while trapping the particles is 



trap 



In 



Znjv) 
Zn{V) 



In 



N 



n\{N 



ujie 



-NE 



and the work extracted when the trapping potentials are turned off is 



off 



In 


\Z{V)] 


= ln 




_Zn{v)_ 





1 

iV! 



V 



V 



N 



-,NE 



Summing (EI]) and ([22 

Wn = Wtrap 



we find the total work to be 
" ]^n\{N -n)\ 



off 



In 



(21) 



(22) 



(23) 



which is independent of E and is equal to the change in uncertainty /„ in ( IT7|) . This 
protocol is optimal and feedback reversible; run in reverse the protocol confines exactly 
n particles in the left half with certainty. 

At first it may be surprising that work can be extracted from this protocol, since 
we are mearly adding and then removing potential minima. However, net work can be 
extracted, since the work extracted while slowly turning on or off a trapping potential 
depends on the total volume accesible to the particles. To see this, consider the simplest 
scenario of turning off one trapping potential with one particle confined to a box of 
volume V . As the depth of the potential minimum becomes shallower, work is done on 
the particle until it escapes from the range of the trapping potential. Once the particle 
leaves, turning off the potential requires no additional work until the particle returns. 
The time for the particle to return depends on the size of the box. For a box of larger 



Designing optimal feedback engines 



12 



volume, the time to return is longer, and the process requires less work. Going back 
to the A^-particle protocol, the work extracted while turning on the trapping potentials 
after the partition has been inserted - when the available volume for each particle is 
V/2 - is more than the work done during the final step as the trapping potentials are 
removed, because the volume V available for the particles to explore is larger. 

When the number of trapping potentials is not equal to the number of particles N, 
this protocol is no longer optimal. The reason being that work can only be extracted 
when a particle can fall into a potential being lowered; the more trapping potentials a 
particle has access to, the more work that can be extracted. If there were less trapping 
potentials then particles, overall less work would be extracted; as there would be fewer 
sites where energy was being removed. If more than trapping potentials are lowered, 
we are able to extract additional work. However, after the partition is removed, each 
particle can explore an even greater number of trapping potentials; the work to turn off 
the potentials would exceed that extracted by turning them on. 

4. Conclusion 

Feedback-reversible processes are optimal, converting all the information acquired 
through feedback into work. In this article, we formulated a strategy, called preparation, 
for designing a collection of optimal protocols given a measured physical observable. In 
the preparation method, optimal protocols are selected by searching for an external 
parameter protocol whose time-reversal prepares the post-measurement state. To 
highlight the utility of the preparation method, we applied it to two pedagogical 
examples - a two- and A^-particle Szilard engine - exhibiting a distinct collection of 
optimal protocols for each. In both examples, we addressed the simplest scenario of 
error-free measurements. When there are measurement errors - for example, if in the 
A^-particle Szilard engine (section 13. 2p . there were a chance to miscount the number 
of particles in the left half of the box - the preparation method still provides a useful 
procedure for selecting an optimal protocol. Furthermore, each of our optimal protocols 
contained at least one infinitely slow step. This is unavoidable as the process must be 
reversible before and after any measurements. Consequently, our method does not 
strictly apply to finite-time processes. However, the preparation method may still 
provide insight into the design of optimal finite-time processes, since an optimal finite- 
time protocol, roughly speaking, is as close to reversible as possible [6l |2^ . 

Generally, we expect the preparation method to be of use whenever the external 
parameter protocol forces a symmetry breaking in the system prior to the measurement, 
such as the insertion of the partition in the Szilard engine. Consider a thermodynamic 
process V during which a system is driven from an initial equilibrium state A through a 
critical point, where the system chooses among several phases or macroscopic states 5, 
with probability pi. In addition, suppose there exists a collection of processes during 
which the symmetry is broken forcibly (not spontaneously), driving the system from A 
to Bi with probability one. Then, according to our recipe this spontaneous symmetry 
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breaking transition can be exploited using the following optimal feedback protocol: start 
in state A, execute process V, measure which state Bi resulted from the symmetry 
breaking, and then run the corresponding process in reverse to drive the system back 
to its initial state A. By construction, this process prepares the post-measurement state 
with unit probability, and therefore extracts as work (W) = —kT J2i Pi logpi, which is kT 
times the information gained in the measurement, (/) = —J2iPi^ogpi. One interesting 
instance of this setup is the Ising model, where a measurement of the system's total 
magnetization after the symmetry breaking phase transition between the paramagnetic 
and ferromagnetic states can be exploited to extract work. This information can be 
utilized by modifying an external magnetic field, as demonstrated in [22]. 

In the introduction, we outlined a general procedure for preparing a collection of 
optimal protocols, original presented in [6l [151 UHl [IZl [HI HH], in which the Hamiltonian 
is instantaneously changed immediately following the measurement in order to make the 
post-measurment state an equilibrium Boltzmann distribution, followed by a reversible 
switching of the external parameters to their final values. These protocols prepare 
the post-measurement states; as such this generic procedure is a special case of the 
preparation method developed here. Though, the implementation of the preparation 
method can lead to a wider variety of protocols. Take for example the two-particle 
Szilard engine discussed in section 13. 1[ Imagine we make a measurement and find 
outcome B, where each particle is confined to a separate half of the box. Let psiz) 
denote the phase space density conditioned on this measurement outcome. In the generic 
procedure, immediately after the measurement we would change the Hamiltonian to 
Hb{z) = — In pb{z), which is a strange Hamiltonian that assigns infinite energy to 
configurations where both particles are in the same half of the box. In contrast, the 
preparation method led to a physically realizable protocol, in which we vary the size of 
the box. 

Finally, we formulated the preparation method only for classical systems. Though, 
the second law of thermodynamics for discrete feedback was originally predicted for 
quantum evolutions [H]. Its mathematical structure resembles the classical version, 
which suggests that feedback-reversible processes are also optimal quantum feedback 
protocols and that the preparation method would also apply to quantum feedback 
engines. Applications of the preparation method to quantum systems holds interesting 
possibilities. For example, in both the classical multi-particle Szilard engines analyzed 
here, the optimal protocols required repulsive particle interactions. In a quantum multi- 
particle Szilard engine composed of fermions, the Pauli exclusion principle induces a 
repulsive interaction of purely quantum origin, which could be exploited to develop a 
collection of optimal feedback protocols. 
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Appendix A. Partition functions for two square hard-core particles in a 
two-dimensional box 

In this appendix, we report the configurational partition functions employed in Sect. 13.11 
for a gas composed of two square particles of width d with hard-core interactions confined 
to a two-dimensional box of width Lx and height Ly. The partition function for hard-core 
particles is the number of distinct microscopic configurations subject to the constraint 
that the centers of the particle be separated by a distance of at least d. In addition, the 
particles have a hard-core interaction with the walls enclosing the box, with the result 
that the center of each particle must be at least a distance d/2 from the edges of the 
box. 

Two partition functions are utlized in our analysis in section 13.11 The first is the 
partition function for the equilibrium state when each particle is confined to separate 
box of dimensions x Ly-. 

rL^-d/2 rLy-d/2 rL^-d/2 rL^-d/2 

Z2{Lx,Ly)= / dxi / dyi / dx2 / dy2 (A.l) 

Jd/2 Jd/2 Jd/2 Jd/2 

= {L,-df{Ly-df. (A.2) 

The second is for the equilibrium state when both particles are confined to the same 
box of dimensions x Ly. This partition function can be expressed as the integral 

1 rL:,~d/2 rLy~d/2 rL^~d/2 rLy~d/2 

Z2{Lx,Ly)= - dxi dyi / dx2 / dy2 

2 Jd/2 Jd/2 Jd/2 Jd/2 

X [e(|xi - X2I -d) + 0(|yi - 2/2I -d)- e(|xi - X2I - d)e{\yi - 2/2I - d)], 

where 9(x) is the Heaviside step function and the preceding factor of 1/2 is included 
because the particles are indistinguishable. The calculation of the above integral can 
be performed using standard methods of integral calculus, with the result, assuming 
Lx > 2d, 

( \{Lx - 2df{Ly - 2df + 2d{Lx - 2d){Ly - 2d){Lx + Ly - Ad) 
Z2iLx, Ly) = I +d^ [{Lx - 2d)2 + {Ly - 2d)'] , Ly > 2d .(A.3) 

[ l{Ly - d)\Lx - 2d)', d<Ly<2d 
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